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DARK-FIELD PHASE CONTRAST IMAGING 

Applicant claims, under 35 U.S.C. § 1 19(e), the benefit of priority of: 1) the 
filing date of December 28, 2000, of U.S. Provisional Patent Application Serial 
Number 60/258,851, filed on the aforementioned date, the entire contents of which are 
5 incorporated herein by reference, and 2) the filing date of February 28, 2001, of U.S. 
Provisional Patent Application Serial Number 60/272,354, filed on the 
aforementioned date, the entire contents of which are incorporated herein by 
reference. 

BACKGROUND OF THE INVENTION 

10 Field of the Invention 

J5 The present invention relates to refractive imaging in general and x-ray 

j jf refractive radiography in particular. 

% Discussion of Related Art 

7 It is well known to use X rays for imaging the internal features of objects in 

1 5 those cases when the object is opaque in the visible optics domain, or when extremely 
fij high spatial resolution is necessary. Traditional x-ray imaging techniques are based 

£ on the absorption contrast, i.e., on the variation of the absorption factor of different 

p " parts of a sample. Therefore, the only way to increase a contrast of small objects in 

traditional x-ray images, is to increase the intensity of an x-ray beam. But this way is 
2 0 unacceptable in many cases. For example, strong x-ray beams cannot be used for 
visualization of the inner structure of integrated circuits because of their possible 
radiation damage, or for medical radiography for safety reasons. The refractive 
contrast, originating from the variation of the refractive indices of different parts of a 
sample, produces far more detailed images of the samples with small features. This 

2 5 type of x-ray imaging is commonly referred to as a phase contrast imaging (PCI). 

However, the direct beam, carrying practically no information about the object, if the 
latter is transparent to x rays, deteriorates the image, bringing additional noise into it. 
Therefore, the direct beam is undesirable. 

Accordingly, it is an object of the present invention to significantly suppress 

3 0 the intensity of a direct beam during phase contrast imaging. 



SUMMARY OF THE INVENTION 

One aspect of the present invention regards an imaging system that includes a 
radiation generator that generates a beam of radiation along a first direction and an 
object that receives the beam of radiation, wherein a first portion of the beam of 
radiation is transmitted through the object along the first direction and a second 
portion of the beam of radiation is refracted along a second direction. An analyzer 
receives the first and second portions of the beam of radiation, the analyzer suppresses 
the intensity of the first portion of the beam of radiation and transmits the second 
portion of the beam of radiation. A detector system that receives from the analyzer 
the suppressed first portion of the beam of radiation and the transmitted second 
portion of the beam of radiation and generates an image of the object. 

A second aspect of the present invention regards a method of imaging an 
object that includes subjecting an object to a beam of radiation that is directed along a 
first direction and analyzing a first portion of the beam of radiation that is transmitted 
through the object along the first direction so that the intensity of the first portion is 
suppressed. Analyzing a second portion of the beam of radiation that is refracted from 
the object. Generating an image of the object based on the suppressed first portion of 
the beam of radiation and the second portion of the beam of radiation. 

An advantage of each aspect of the present invention is the reduction of the 
negative influence of the direct x-ray beam in an x-ray refractive imaging system. 

A second advantage of each aspect of the present invention is the reduction in 
radiation flux that an object is subjected to during imaging. 

Additional objects and advantages of the invention will become apparent from 
the following description and the appended claims when considered in conjunction 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 schematically shows the fundamentals of x-ray refractive radiography 
in general; 

FIG. 2 schematically shows a top view of a first embodiment of an imaging 
system according to the present invention; 



FIG. 3 schematically shows a side view of the imaging system of FIG. 2; 

FIG. 4 schematically shows a top view of a second embodiment of an imaging 
system according to the present invention; 

FIG. 5 schematically shows a side view of the imaging system of FIG. 4; 

FIG. 6 schematically shows a mode of scanning to be used with the imaging 
system of FIG. 4; 

FIG. 7 schematically shows a side cross-sectional view of a first embodiment 
of an analyzer to be used with the imaging systems of FIGS. 2 and 4 according to the 
present invention; 

FIG. 8a shows angular reflection and absorption curves for the analyzer of 

FIG. 7; 

FIG. 8b shows an enlarged portion of the reflection curve of FIG. 8a; 

FIG. 9 schematically shows the intensity distribution of the analyzer of FIG. 7; 

FIG. 10 schematically shows a test object that can be imaged by the imaging 
systems of FIGS. 2 and 4; 

FIGS. 1 la-c show pictures of the test object of FIG. 10 when imaged by the 
imaging systems of FIGS. 2 and 4 but without the analyzer of FIG. 7; and 

FIGS. 12a-c show pictures of the test object of FIG. 10 when imaged by the 
imaging systems of FIGS. 2 and 4 and with the analyzer of FIG. 7. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The general principles of x-ray refractive radiography are shown in FIG. 1. In 
particular, a parallel x-ray beam 100 with a wave vector k 0 penetrates the object 102, 
transparent to x-rays. Because of the refraction of x-rays on the inner structure of the 
object 102, the output beam 104 is composed of the original wave 106 with the vector 
k 0 and the refracted waves 108 with vectors k(9) slightly declined from the initial 
vector k 0 . An analyzer 1 10 then differentiates the output waves 104 by their direction 
9 . Thus, the resultant image registered in the image plane 112 directly behind the 
analyzer 110 is formed by the intensity variations R(6) proportional to the angle of 
refraction 0 . Therefore, the resultant intensity distribution in the image plane 112 
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may be presented in the form of a sum 

I(r) = R[9(r)] + Io (1) 
where I 0 is the intensity of the direct beam 106, which may be much greater than the 
useful term R[9(r )]. The direct beam 106 carries practically no information about the 
5 object, if the latter is transparent to x rays, and at the same time it deteriorates the 
image, bringing additional noise into it. Therefore, the detector shot noise produced 
by the direct beam decreases the signal-to-noise ratio, which is of a primary 
importance when an image is recorded electronically. For these reasons, the direct 
beam 1 06 is undesirable. 
10 ° ne embodiment of an imaging system that reduces the effects of the direct 

jz beam 106 is shown in FIGS. 2 and 3. In this embodiment, an x-ray tube 114 that 

Q works in the linear projection mode is used. As shown in FIG. 2, the x-ray tube 1 14 

generates a beam 1 1 5 so that the long side of the focus 1 1 6 of the beam 1 1 5 is in the 
w plane of incidence. The beam 115 is directed to a monochromator 118 that may be 

in 15 composed of two crystals 120, 122 that are well known in the art. The two crystals 
^ 120, 122 are selected so that they strongly disperse the beam 115 so as to generate 

ill highly parallel x-ray beams 1 00. The width of the beam 1 00 in the plane of incidence 

m is limited by the length of either the x-ray tube focus 1 1 6 or the mirror 1 1 0 while the 

rt * 

g height of beam 100 is determined by the vertical divergence Y of the beam after the 

2 0 monochromator 1 1 0 according to the equation: 



y « 
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Here a and Act are the angle of incidence and the width of the Bragg reflection for 
the monochromator crystals 120, 122. The value of J typically does not exceed 3°. 

In the embodiment of FIGS. 2 and 3, the object 102 is preferably no larger than 
2 5 several millimeters so that the object 102 is fully covered by the x-ray beam 100. 
Accordingly, there is no need to move the object 102 during imaging. Please note that 
the divergence of the probe beam 100 should not exceed approximately 10 arc 
seconds. Apart from the divergence, the spectral width of the probe beam 100 also 
plays a significant role because the mirror 110's resonant angle depends on the 



wavelength. In the case of the x-ray tube 1 14 emitting X-ray beams about the CuKa 
doublet, the natural relative spectral width of the CuKa doublet including side wings 
is of the order of 4x1 0~ 3 . Such spectral widening increases the angular width of the 
resonance by a value of approximately 10 arc seconds. Therefore, the probe beam 100 
preferably should not only have an extremely small divergence, but good spectral 
purity as well. Coupling the beam 115 with a double crystal monochromator 118 
working in a strongly dispersive mode with both Si (111) reflections achieves such 
divergence and spectral purity. Using such a double crystal monochromator, the 
horizontal divergence of the probe beam 100 is approximately equal to 7.4 arc 
seconds, while the relative spectral width is approximately 5xl0" 4 , i.e., an order of 
magnitude less than the natural spectral width of the CuK^ radiation. Thus, the effect 
of spectral widening can be neglected. 

After the beam 100 interacts with the object 102, the beam 104 is directed to 
an analyzer 110 that suppresses the intensity of the original wave or beam 106 by 
several orders of magnitude in a manner as schematically shown in FIG. 9. The 
suppressed beam 106 and the refracted beam 108 are directed to the imaging plane 
112 where a detector, such as an x-ray charge coupling device (CCD) 113, receives 
the beams. The detector then sends a signal to a processor (not shown) that generates 
an image that is formed on a display (not shown). 

One embodiment of an analyzer 110 that can suppress the intensity of the 
beam 106 is shown in FIG. 7. In particular, the analyzer 110 of FIG. 7 is a specially 
designed multilayer mirror 124. The reflective coating of the x-ray multilayer mirror 
124 is composed of many altering layers of materials with large and small atomic 
numbers. For instance, the layers 126 with large atomic numbers may be made of 
tungsten while the layers 128 with small atomic numbers may be made of boron- 
carbide, i.e., B 4 C. The thickness of the layers may differ, but they are typically of the 
order of 10A-50A. The interfacial roughness is equal to 5A. 

As described in "X-Ray Multilayer Mirrors with an Extended Angular Range," 
by Protopopov et ah, Optics Communications Vol. 158 (1998), pp. 127-140, the entire 
contents of which are incorporated herein by reference, it is possible to control the 
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shape of the angular and spectral reflection curves by altering the thickness of the 
layers 126 and 128. Depth-graded multilayer mirror can be analyzed as a combination 
of traditional multilayer structures with constant periods, imposed one on another and 
calculated for different grazing angles. The range of periods is chosen so as to cover a 
5 given angular region. At a definite grazing angle only one structure reflects, namely 
that whose period matches the grazing angle. Others, on top are transparent to x-rays, 
although introducing some attenuation. Varying slightly the thickness of layers it is 
possible to make the partial reflected waves approximately counterphased at a specific 
grazing angle 9 , so as to obtain as small reflection at this angle as possible. 
10 Moreover, the total reflection can be made even less if not only the phases of the 
p partial waves are opposite each to another, but the coming and reflected waves 

H produce interference pattern whose maxima at this particular angle coincide with the 

yy 

Ul layers of heavy material, introducing additional absorption. Thus, it is possible to 

ffj design a mirror with deep (the reflectivity of the order of 10" 2 -10" 3 ) and narrow 

- Pi 

15 (several arc seconds) resonant gap in the angular reflection curve as shown in FIGS. 
8a-b. The roles of reflection and absorption are clear from the solid and dashed 
fy curves, respectively, in FIG. 8a. In addition, the sensitivity of the scheme with respect 

5; to the refracted beams 108 is determined by the sharpness of the reflection curve 

around the resonant angle 0 r . 
2 0 It should be kept in mind that the angular divergence co of the beam 100 is a 

very important parameter determining both the depth and the angular width of the 
resonant gap. That is clear from FIG. 8b, where three reflection curves labeled 1, 2, 3, 
calculated for different values of co (curve 1 : 2a>=0; curve 2: 2co=0.003° (10 arc sec); 
curve 3:- 2a>=0.016° (1 arc min.)), are shown. 

2 5 As shown in FIGS. 12a-c, the sharpness of the gap in the reflection curve of 

the multilayer mirror 124 is sufficient to effectively detect small-contrast images. In 
particular, when a test object 102 is in the form of a carbon sphere of 0.8mm in 
diameter, surrounded by a layer of water 1mm thick as shown in FIG. 10 and is used 
as the object of the imaging system 101 of FIG. 2, then improved images are obtained. 

3 0 The three images of FIGS. 12a-c are obtained in a geometrical optics approximation, 



for CuKa radiation (1.54A), using the reflection curve 2 in FIG. 8b. The divergence of 
the primary beam 100 is chosen to be 2co=0.003° (10 arc sec), that is a practically 
attainable value, while the distance between the object 102 and the image plane is 
100mm. Under these conditions the loss of resolution due to the beam divergence is 
negligible. The images of FIGS. 12a-c differ from each another by the total number of 
photons/flux in the primary beam 100 where FIGS. 12a-c have fluxes of 10 6 
photons/mm 2 , 10 5 photons/mm 2 and 10 4 photons/mm 2 , respectively. In addition, the 
horizontal axes of the images lie in the incidence plane of the multilayer mirror 124. 
Therefore, the contours of the sphere are disconnected at the top and at the bottom, 
where the projections of the refracted beam vectors onto the incidence plane have no 
deflections. 

As shown in FIGS. 1 la-c, the test object of FIG. 10 was subjected to x-rays 
under the same conditions as FIGS. 12a-c, respectively, so as to form 
absorption-contrast images simulated for a traditional x-ray radiography scheme, i.e., 
without the analyzer. Comparing the images shows that contour of the carbon sphere 
pattern is recognized in each of the images of FIGS. 1 la and 12a obtained with a 
photon flux of 10 6 photons/mm 2 . However, at lower fluxes the contours and the 
pattern itself can hardly be recognized in absorption-contrast images (FIGS. 1 lb, 
11c). As far as the refraction-contrast images of FIGS. 12b and 12c are concerned, 
the sphere pattern can be recognized even at a flux of 10 4 photons/mm 2 . Thus, the 
radiation dose can be decreased by approximately two orders of magnitude in 
comparison with traditional x-ray radiography. This is advantageous when imaging 
biological objects. 

If it is desired to image objects that are larger than 2mm and have dimensions 
up to 150-200 mm, then a modified imaging system can be employed. This is 
advantageous for biological and medical applications. An embodiment of such an 
imaging system is shown in FIGS. 4-6. In this embodiment, the x-ray tube 1 14 works 
in the point projection mode. The width of the beam in the plane of incidence is 
limited by the x-ray tube focus, and is an order of magnitude less than that for the 
imaging system of FIGS. 2-3. Consequently, the length of the mirror 1 10 in this 
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direction may be much less than in the previous case. 

As shown in FIG. 4, the x-ray tube 114 generates a beam 115 that is directed to 
the monochromator 118 that is composed of two crystals 120, 122 that are similar to 
those described previously with respect to the imaging system of FIG. 2. Again, the 
two crystals 120, 122 are selected so that they strongly disperse the beam 1 15 so as to 
generate highly parallel x-ray beams 100. 

In the embodiment of FIGS. 4 and 5, the object 102 is preferably larger than 
the width of the x-ray beam 100. Accordingly, there is a need to move the object 102 
relative to the detector 113 during imaging as shown in FIG. 6. The object 102 may 
be moved by a number of conventional devices, such as stepping motor driven 
translation stages that are well known in optical imaging. 

After the beam 100 interacts with the object 102, the beam 104 is directed to 
an analyzer 110 that suppresses the intensity of the original wave or beam 106 by 
several orders of magnitude in a manner as schematically shown in FIG. 9. The 
suppressed beam 106 and the refracted beam 108 are directed to the imaging plane 
112 where a detector, such as an x-ray charge coupling device 113, receives the 
beams. The detector then sends a signal to a processor (not shown) that generates an 
image that is formed on a display (not shown). The analyzer 110 preferably has a 
structure that is similar to that as the analyzer 110 used in the imaging system of 
FIGS. 2-3. 

As shown in FIG. 6, the object 102 is scanned in the plane of incidence in the 
direction transversal to the x-ray beam 100, so that each moment of time only a small 
fraction of the object is investigated. During each moment of time t the detector 
signal can be described by the matrix tyy(f), where ' and j are the ordinal numbers 
of its sensitive elements. The signals corresponding to the same row j but different 
column / differ each from another only by the time delay equal to /T , where X i s 
the time interval during which the object is shifted by a distance equal to a single 
detector element. Therefore, it is possible to average the signals from different 
columns if only take into account the delay. Such an averaging will rise the sensitivity 



and signal-to-noise ratio because the noise in the channels is uncorrected. The time t 
corresponds to the first discrete coordinate of the image k by the formula t = kx 9 
while the second discrete coordinate of the image is the row number j . Thus, the 
averaged discrete image can be written in the following form: 

f "/=l fn i=l 
where m is the number of columns in the detector. 

The images obtained using the imaging system of FIGS. 4-6 are comparable to those 
images obtained using the imaging system of FIGS. 2-3 and thus the imaging system 
of FIGS. 4-6 enjoys many of the advantages of the imaging system of FIGS. 2-3 as 
described previously. 

While the above description constitutes the preferred embodiments of the 
present invention, it will be appreciated that the invention is susceptible of 
modification, variation and change without departing from the proper scope and fair 
meaning of the accompanying claims. 



